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We consider an electronic spin, such as a nitrogen-vacancy (NV) center in diamond, weakly coupled 
to a large number (bath) of nuclear spins, and subjected to the Rabi driving with a periodically 
alternating phase (multiple rotary echo). We show that by switching the driving phase synchronously 
with the precession of a given nuclear spin, the interaction to this spin is selectively enhanced, while 
the rest of the bath remains decoupled. The enhancement is of resonant character. The key feature 
of the suggested scheme is that the width of the resonance is adjustable, and can be greatly decreased 
by increasing the driving strength. Thus, the resonance can be significantly narrowed, by a factor of 
10-100 in comparison with the existing detection methods. Significant improvement in selectivity is 
explained analytically and confirmed by direct numerical many-spin simulations. The method can 
be applied to a wide range of solid-state systems. 

PACS numbers: 76.70.-r,61.05.Qr,76.30.Mi,03.67.Lx 


Detecting a single nuclear spin is an ultimate goal in 
nuclear magnetic resonance, actively pursued using var¬ 
ious techniques EH- Moreover, one-by-one detection, 
characterization, and manipulation of individual nuclear 
spins in solids is vital for harnessing them as a resource 
for quantum information processing fol-|20i|. The crucial 
problem, hindering the progress in this area, is to sepa¬ 
rate the weak signal, produced by a given weakly coupled 
nuclear spin, from the much stronger background created 
by all other spins and by the ambient magnetic noise. 
I.e., the detector must possess high selectivity to sepa¬ 
rate the signal of the target spin from the action of other 
spins and from magnetic noise. Thus, the focus of the 
present work is to develop a method for highly selective 
and sensitive one-by-one detection and manipulation of 
individual nuclear spins in solids. 


It has been recently shown that an electron spin of the 
nitrogen-vacancy (NV) centers in diamond is an excel¬ 
lent nuclear spin detector 21-271. By applying the 


dynamical decoupling protocols to the NV spin (pulse 
[5|-l7| or continuous [S), 127]), the individual nuclear spins 
around the NV center can be detected one by one in a 
resonant manner: the NV spin becomes decoupled from 
all nuclear spins, except the one which satisfies a strin¬ 
gent resonance condition. The resonating nuclear spin 
strongly affects the NV spin motion, and hence can be 
detected, and the width of the resonances produced by 
different nuclear spins determines selectivity. However, 
in all current schemes this width is naturally limited, be¬ 
ing determined by the strength of the hyperfine coupling 
between the target nuclear spin and the NV spin. 


Overcoming this limitation, and drastically narrowing 
the resonances in a systematic way, in order to signifi¬ 
cantly improve the detection selectivity, is of much im¬ 
portance. Here we present a scheme which achieves that 
goal. It employs the periodically changing Rabi driv¬ 


ing on the NV center’s spin (multiple rotary echo) [2g 
l30j with a specially chosen period to detect the nuclear 
spins. The width of a resonance, corresponding to detec¬ 
tion of a given nuclear spin is inversely proportional to 
the Rabi driving strength, and therefore can be directly 
adjusted in experiment. As a result, the resonance peak 
width is narrowed by a factor of 10 — 100 in compari¬ 
son with the existing schemes. This narrowing greatly 
enhances the resolution, which is very important for ac¬ 
curate characterization of the nuclear-spin environment 

, and for NV-based nuclear 


of the NV center 


magnetic resonance at nanoscale 21, mm. Thus 


our 

approach combines excellent protection from magnetic 
noise offered by the rotary echo, and the greatly enhanced 
selectivity in detection of individual nuclear spins. This 
scheme can also be used with other spins as detectors, 
and in some conventional ensemble magnetic resonance 
experiments involving electronic and nuclear spins. 


We consider the electronic spin S' = 1 of the NV cen¬ 
ter, subjected to a moderate (tens to hundreds of Gauss) 
static bias field along the symmetry axis. The NV spin 
possesses three well-separated states uinv = 0,1 and 
— 1 (denoted below as |0,±1}), and is manipulated by 
applying the microwave driving at resonance with the 
transition between the states |0) and |1) (Fig. QJi,); the 
level | — 1) remains idle and will be further ignored. The 
NV spin is coupled to a nuclear spin I: e.g., a 13 C nu¬ 
clear spins in diamond (I = 1/2), coupled to the NV 
center via dipolar interaction. To simplify consideration, 
we temporarily omit the coupling of the NV spin to the 
other sp ins, considering many nuclear spins below (see 
also [311 ]). In the coordinate frame rotating with the fre¬ 
quency of the resonant microwave driving, the system is 
described by the Hamiltonian 


H — hS x + (A||/ z -f- A±I X )S Z + [lol + A||/2)/ z , (1) 
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where S x = 5 (| 1)<0| + |0)(1|) and S z = 5 (|1)<1| - |0)(0|) 
describe the NV spin, and I x>z are the operators of the 
target nuclear spin, h is the Rabi driving field (usually, 
few MHz to few tens of MHz), u>l is the Larmor frequency 
of the nuclear spin, Ay = Uh cos 6 and A± = uu sin 6 are 
the components of the hyperfine coupling |31| , and we 
assume h = 1 throughout the text. 

During the rotary echo experiment, the direction of the 
Rabi driving field (i.e. the phase of driving) periodically 
changes between +x (along the £-axis) and —x (oppo¬ 
site to the x-axis), as shown in Fig.[lji. Here we consider 
symmetrized protocol, consisting of N rotary echo cycles; 
within each cycle the driving field is along +£ during the 
two outer segments, each of duration T, and is along —x 
during the inner segment of duration 2 T. Therefore, the 
Hamiltonian m describes the two outer segments of the 
rotary echo cycle, and the Hamiltonian for the inner seg¬ 
ment is obtained by replacing h —> —h. Below, we denote 
these two Hamiltonians as H + and ff_, respectively, so 
the evolution operator for a single cycle is 


U = exp (~iTH+) exp(-2 iTH_) exp (~iTH+), (2) 


and the full iV-cycle protocol is described by the evolution 
operator U(N) = XJ N . Generally, periodic switching of 
the Rabi field in the regime of strong driving (h 0Jh) 
leads to highly efficient decoupling of the NV electronic 
spin from the surrounding spins [ 2 bl. fiif)L Ebj| . However, in 
a special resonant regime, where the Rabi field is switched 
synchronously with the Larmor precession of the target 
nuclear spin, quantum states of the NV center spin and 
the target nuclear spin become entangled. After many 
rotary-echo cycles, N 1, this entanglement becomes 
detectable despite the weak coupling between the NV 
and target spins. 

To qualitatively understand this phenomenon, we an¬ 
alyze the evolution of the system assuming small cou¬ 
pling (ujh -C h), so that the zero-order Hamiltonians are 
H±^ = ±hS x + (ul + A\\/2)I Z . The corresponding zero- 
order evolution operator is given by 

c/(°) = exp[— i4T(iOL + A||/ 2 )J Z ] = cos 2 ip — 2 il z sin 2 ip, 

( 3 ) 

where <p = (ojl + A\\/2)T. I.e., the nucleus and the NV 
center are practically decoupled: over one cycle, the nu¬ 
clear spin rotates around the z-axis by the angle 2 ip, and 
this rotation does not depend on the state of the NV spin. 
However, if the driving is switched in resonance with the 
nuclear spin precession, so that the switching time has a 
special value, 


7r(2fc — 1) 

2 wl ± A|| 


k = 1 , 2 ,.., 


( 4 ) 


then sin 2ip = 0 , and the operator U ^ is unity, i.e. the 
effect of the zero-order Hamiltonian is null after one cy¬ 
cle. Then the smaller higher-order corrections will be¬ 
come important, since their effect will accumulate over 


m NV~ 1 ,-\ 
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FIG. 1: (Color online) (a) Rotary echo protocol. Rabi driv¬ 
ing is applied to the NV electron spin at the frequency of the 
transition between the states m,NV = 0 and ulnv = 1. The 
phase of driving is periodically switched by 180°, i.e. the driv¬ 
ing field in the rotating frame periodically changes between 
+h and —h. (b) Generally, the rotary echo protocol efficiently 
decouples the nuclear spin from the NV center: evolution of 
the nuclear spin is practically independent of the NV spin 
state (left). However, when the switching is on resonance 
with the nuclear spin precession {T = T*,), the nuclear spin 
entangles with the NV spin: its evolution is conditioned on 
the initial NV spin state (right), (c) When T becomes equal 
to the resonance value Tk, entanglement with the nuclear spin 
leads to faster decay of the NV electron spin, seen as a reso¬ 
nant drop of the signal. Graph shows the rotary echo signal 
(defined everywhere in the paper as 2(S Z (N))) after N = 50 
cycles as a function of T, for lol = 2-ir ■ 430 kHz, driving 
h = 2-7T-3.4 MHz, and coupling constants Ay = 27T-—9.55 kHz, 
A± = 2-k ■ 38.2 kHz. The resonances of the orders k = 1 and 
k = 2 are marked with green dashed lines. 


many rotary echo cycles, undisturbed by the zero-order 
terms. The detailed analysis [31] shows that in this reso¬ 
nance case the nuclear spin evolution changes drastically: 
the rotation axis tilts towards y axis, the rotation angle 
per cycle is of order of A±/h , and most importantly, the 
nuclear spin entangles with the NV spin through the op¬ 
erator S y = |(|0)(1| — |1)(0|). I.e., at resonance the 
nuclear spin rotates in different directions depending on 
whether the NV spin is in the state |+) or |—), where |±) 
are the eigenstates of S y with the eigenstates s y = ±1/2, 
respectively. 

We consider an experiment where the NV spin is pre¬ 
pared in the state uinv = 1, and the operator S Z (N) is 
measured after N rotary echo cycles. At resonance, the 
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initial coherent superposition |1) = (|+) + |—))/\/2 is de¬ 
cohered in the basis |±) as a result of entanglement with 
the unpolarized nuclear spin, so that the resonance can 
be detected by sudden onset of strong decay of S Z (N) 
Analytical solution, valid for long times (N 1), is de¬ 
rived in Supplementary Material [31 ]: 


Signal = 2{S Z (N)) = 1 — , , 2 9 L(N), (5) 

and in the vicinity of the resonance g/J. ~ 
( A±/h) sin p cos 4> and L(N) « sin 2 (N gfT), where 
<j> = hT. Thus, the signal as a function of the switch¬ 
ing time T (or the angle ip) has Lorentzian shape 
with the width Aip = (A±/h) cos (f> and the depth 
2 sin 2 [N(A±/h) cos </>]. To maximize the depth, the 
driving should be adjusted to have </> = hT = wm 
with integer m. This also minimizes interference be¬ 
tween different nuclear spins and the influence of the 
fluctuations in the driving power, see Supplemental 
Information & and below we always assume this 
condition satisfied. Typical dependence of the signal 
on T is shown in Fig. [TJ;, where the resonances at Tk- i 
and Tfc =2 are clearly seen. The analytical solution ([5]) 
is very precise: it practically coincides with the direct 
numerical simulation of the two-spin system described 
by the Hamiltonian ©■ 

For a weakly coupled nucleus, located sufficiently far 
from the NV center, the hyperfine coupling is deter¬ 
mined primarily by dipolar interactions, so that A\\ = 
( 7 e 7 n/^? 3 )[l — 3 Rl/R 2 }, where R is the vector connect¬ 
ing the positions of the NV center and the nucleus, and 
7 e and are the electronic and nuclear gyromagnetic 
ratios, respectively. Therefore, the nuclei from different 
locations produce the peaks at different values of T, and 
hence can be resolved. 

To test our approach, and to illustrate its performance 
under realistic circumstances, we performed direct sim¬ 
ulations of the rotary-echo detection for the NV center 
coupled to 14 nuclear spins of 13 C randomly located in 
diamond lattice. The simulation results, given in Fig. [2] 
for k = 3, clearly show the sharp, well-resolved peaks, 
corresponding to different nuclei, with selectivity in sub- 
kHz region, even for very modest driving h ~ 27r ■ 5 MHz 
(see also Supplemental Materials for more details). Some 
nuclei, located at the symmetry-related positions in the 
lattice, have the same A^ and hence the same resonance 
value of T; to resolve them, the static bias field should 
be tilted away from the symmetry axis of the NV cen¬ 
ter jb]. The simulations also confirm our conclusion (31] 
that the NV-mediated interaction between the nuclear 
spins, caused by the driving, does not noticeably affect 
the detection efficiency. 

For a fixed total interrogation time tjv = 4 NT, the res¬ 
onance width determines selectivity [32j of the protocol, 
i.e. the ability to resolve two nuclei with close hyperfine 
couplings and u b h (i.e. ~ u> h , |- u b h \ <C uj h ). 
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FIG. 2: (Color online) (a) Rotary echo signal after N = 100 
cycles, as a function of the switching time T for a NV spin cou¬ 
pled to 14 13 C nuclear spins, randomly placed in the diamond 
lattice with natural abundance, for ujl = 2n ■ 428 kHz. Driv¬ 
ing h is adjusted so that hT = 287r (so that 4.59 < (h/2n) < 
5.28 MHz). The region corresponding to the resonances of the 
order k = 3 is shown. Theoretically expected resonances are 
marked with green lines; each of the lines marked with arrows 
correspond to a pair of nuclei located at the symmetry-related 
sites, and having almost the same value of Ay. Parameters of 
the nuclear spins are given in Ref. l3ll . (b)-(e) Magnified view 
of the different areas of the panel (a). 


High selectivity of the suggested scheme is determined by 
the small width of the resonances, AT = A±/(hu>i,) ~ 
Wft,/(/iWi); here and below we neglect the difference be¬ 
tween A|| and A±, taking A\\ ~ A± ~ w/,. The res¬ 
onance width AT should be compared to the distance 
between the fc-order resonances from two nuclear spins, 
which is ~ (2k — 1 )(w£ — w]()/w|. Thus, the resonances 
are resolved when (wj) — ufy/ujh ~ [ujL/h]/(2k — 1). 
For a typical experiment, detecting 13 C spins, with the 
bias field of 400 G (ujl = 2-7T ■ 428 kHz) and driving of 
h = 27t • 10 MHz, for k = 3, the condition of resolved 
resonances is (oj% — ufy/uh > 1/100. 

The fact that the resonance width can be adjusted by 
simply changing the driving h is a key feature of the 
proposed protocol, which leads to great improvement in 
selectivity. In the currently existing detection schemes 
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|32| . the width of the resonance is naturally limited, be¬ 
ing determined by the coupling constants A\\ and A±. 
The scheme proposed here is the first, to our knowledge, 
where the width of the resonances can be tuned at exper¬ 
imentalist’s will; this allows narrowing of the resonances 
by a factor of 10-100 in comparison with other existing 
schemes (see [31] for details), albeit at the expense of the 
longer interrogation time. Also, stronger driving greatly 
improves the coherence time of the NV spin 29, 30, 33 ]. 
Moreover, using the rotary echo protocol with the opti¬ 
mal choice of driving h significantly reduces the impact 
of slow fluctuations in the driving strength, which is an 
important limiting factor for the Hartmann-Hahn double 
resonance detection [8]. 


Above we omitted the NV center’s own 14 N (or 15 N) 
nuclear spin: we assume that it is 100% polarized in the 
state with a given m/, and that the driving is applied 
at the frequency of the corresponding transition. With¬ 
out such polarization, the hyperfine levels which are not 
in resonance with the Rabi driving field are quickly de¬ 
cohered, since the on-site hyperfine coupling (2.15 MHz 
for 14 N, and 3.03 MHz for 15 N) is comparable to a typi¬ 
cal driving strength, so that overall detection quality de¬ 
grades [31(. Fortunately, polarization of the NV’s own 
nuclear spin can be easily achieved in different ways 
[l2l . ll5ll34H36l ]. and incorporated in the experimental pro¬ 
tocol. 


Concluding, we presented a scheme for one-by-one de¬ 
tection of the nuclear spins weakly coupled to an elec¬ 
tron spin (e.g. the NV center in diamond), using the ro¬ 
tary echo protocol. The key feature of this scheme is the 
ability to experimentally adjust, and thus greatly nar¬ 
row, the width of the resonance peaks corresponding to 
the detected nuclear spins. When compared with the 
existing approaches, the scheme proposed here provides 
significant improvement in the selectivity of the nuclear 
spin detection, albeit being more demanding to the qual¬ 
ity of experimental setup (timing precision, fast switch¬ 
ing of the driving phase, increased interrogation time, 
etc.). Moreover, the scheme discussed here, and the con¬ 
siderations above, are not specific to the NV center, and 
can be applied to other electron spins, or to conventional 
ensemble double electron-nuclear resonance experiments 
[3l| . The main requirements are: (a) well-defined static 
hyperfine couplings with anisotropy component, (b) nu¬ 
clear Larmor frequency sufficiently large in comparison 
to hyperfine couplings, and (c) sufficiently long coherence 
times of the electron and nuclear spins; a wide range of 
solid-state systems satisfy these conditions. 

We thank R. Hanson and T. H. Taminiau, and M. 
Raikh for helpful and important discussions. This work 
was supported by the Department of Energy Basic En¬ 
ergy Sciences under Contract No. DE-AC02-07CH11358. 


* Electronic address: slava@ameslab.gov 

[1] C. L. Degen, M. Poggio, H. J. Mamin, C. T. Rettner, 
and D. Rugar, Proc. Nat. Acad. Sci. 106, 1313 (2009). 

[2] R. Vincent, S. Klyatskaya, M. Ruben, W. Wernsdorfer, 
F. Balestro, Nature 488, 357 (2012). 

[3] J. J. Pla, K. Y. Tan, J. P. Dehollain, W. H. Lim, J. J. 

L. Morton, F. A. Zwanenburg, D. N. Jamieson, A. S. 
Dzurak, and A. Morello, Nature 496, 334 (2013). 

[4] J. J. L. Morton, D. R. McCamey, M. A. Eriksson, S. A. 
Lyon, Nature 479, 345 (2011). 

[5] T. H. Taminiau, J. J. T. Wagenaar, T. van der Sar, F. 
Jelezko, V. V. Dobrovitski, and R. Hanson, Phys. Rev. 
Lett. 109, 137602 (2012). 

[6] S. Kolkowitz, Q. P. Unterreithmeier, S. D. Bennett, and 

M. D. Lukin, Phys. Rev. Lett. 109, 137601 (2012). 

[7] N. Zhao, J. Honert, B. Schmid, M. Klas, J. Isoya, M. 
Markham, D. Twitchen, F. Jelezko, R.-B. Liu, H. Fedder, 
and J. Wrachtrup, Nature Nanotech. 7, 657 (2012). 

[8] P. London, J. Scheuer, J.-M. Cai, I. Schwarz, A. Ret- 
zker, M. B. Plenio, M. Katagiri, T. Teraji, S. Koizumi, 
J. Isoya, R. Fischer, L. P. McGuinness, B. Naydenov, and 
F. Jelezko, Phys. Rev. Lett. Ill, 067601 (2013). 

[9] M. V. G. Dutt, L. Childress, L. Jiang, E. Togan, J. Maze, 
F. Jelezko, A. S. Zibrov, P. R. Hemmer, and M. D. Lukin, 
Science 316, 1312 (2007). 

[10] P. Cappellaro, L. Jiang, J. S. Hodges, M. D. Lukin, Phys. 
Rev. Lett. 102, 210502 (2009) 

[11] W. Pfaff, B. Hensen, H. Bernien, S. B. van Dam, M. S. 
Blok, T. H. Taminiau, M. J. Tiggelman, R. N., Schouten, 
M. Markham, D. J. Twitchen, R. Hanson, Science (2014), 
DOI: 10.1126/science. 1253512 

[12] T. H. Taminiau, J. Cramer, T. van der Sar, V. V. Do¬ 
brovitski, R. Hanson, Nature Nano. 9, 171 (2014). 

[13] P. Neumann, N. Mizuochi, F. Rempp, P. Hemmer, H. 
Watanabe, S. Yamasaki, V. Jacques, T. Gaebel, F. 
Jelezko, and J. Wrachtrup, Science 320, 1326 (2008). 

[14] D. D. Awschalom, L. C. Bassett, A. S. Dzurak, E. L. Hu, 
J. R. Petta, Science 339, 1174 (2013). 

[15] G. D. Fuchs, G. Burkard, P. V. Klimov, D. D. 
Awschalom, Nature Phys. 7, 789 (2011). 

[16] P. C. Maurer, G. Kucsko, C. Latta, L. Jiang, N. Y. Yao, 
S. D. Bennett, F. Pastawski, D. Hunger, N. Chisholm, 
M. Markham, D. J. Twitchen, J. I. Cirac, M. D. Lukin, 
Science 336, 1283 (2012). 

[17] Y. Zhang, C. A. Ryan, R. Laflamme, J. Baugh, Phys. 
Rev. Lett. 107, 170503 (2011). 

[18] J. S. Hodges, J. C. Yang, C. Ramanathan, D. G. Cory, 
Phys. Rev. A 78, 010303 (2008). 

[19] L. Childress, M. V. Gurudev Dutt, J. M. Taylor, A. S. 
Zibrov, F. Jelezko, J. Wrachtrup, P. R. Hemmer, and M. 
D. Lukin, Science 314, 281 (2006). 

[20] V. V. Dobrovitski, G. D. Fuchs, A. L. Falk, C. Santori, 
and D. D. Awschalom, Ann. Rev. Cond. Matter Phys. 4, 
23 (2013). 

[21] R. Schirhagl, K. Chang, M. Loretz, C. L. Degen, Ann. 
Rev. Phys. Chem. 65, 83 (2014). 

[22] A. Ajoy, U. Bissbort, M. D. Lukin, R. Walsworth, and P. 
Cappellaro, arXiv: 1407.3134. 

[23] M. Kost, J. Cai, and M. B. Plenio, arXiv:1407.6262. 

[24] N. Zhao, J.-L. Hu, S.-W. Ho, J. T. K. Wan, and R. B. 
Liu, Nat. Nano 6, 242 (2011). 







5 


[25] H. J. Mamin, M. Kim, M. H. Sherwood, C. T. Rettner, 
K. Ohno, D. D. Awschalom, and D. Rugar, Science 339, 
557 (2013). 

[26] T. Staudacher, F. Shi, S. Pezzagna, J. Meijer, J. Du, C. 
A. Meriles, F. Reinhard, and J. Wrachtrup, Science 339, 
561 (2013). 

[27] J. Cai, F. Jelezko, M. B. Plenio, and A. Retzker, New J. 
Phys. 15, 013020 (2013). 

[28] I. Solomon, Phys. Rev. Lett. 2, 301 (1959). 

[29] A. Laraoui and C. A. Meriles, Phys. Rev. B 84, 161403 

( 2011 ). 

[30] C. D. Aiello, M. Hirose, and P. Cappellaro, Nature 
Comm. 4, 1419 (2013). 

[31] Supplemental Materials 

[32] The protocols discussed in this work are, in fact, two- 
dimensional, with the signal depending on some ad¬ 
justable parameter (T for our scheme, r for the pulse 
dynamical decoupling, etc.), and on the total interroga¬ 
tion time. Ref. [S}| proposes to analyze the latter depen¬ 


dence as a way to increase resolution; this option can also 
be used within the pulse dynamical decoupling detection, 
and within the protocol presented here. However, the re¬ 
sulting improvement strongly depends on experimental 
details, and is beyond the scope of the present paper. 

[33] V. V. Mkhitaryan and V. V. Dobrovitski, Phys. Rev. B 
89, 224402 (2014). 

[34] V. Jacques, P. Neumann, J. Beck, M. Markham, D. 
Twitchen, J. Meijer, F. Kaiser, G. Balasubramanian, F. 
Jelezko, and J. Wrachtrup, Phys. Rev. Lett. 102, 057403 
(2009). 

[35] P. Neumann, J. Beck, M. Steiner, F. Rempp, H. Fedder, 
P. R. Hemmer, J. Wrachtrup, F. Jelezko, Science 329, 
542 (2010). 

[36] T. van der Sar, Z. H. Wang, M. S. Blok, H. Bernien, T. H. 
Taminiau, D. M. Toyli, D. A. Lidar, D. D. Awschalom, R. 
Hanson, and V. V. Dobrovitski, Nature 484, 82 (2012). 


